Neutrino-nucleus cross section in the impulse approximation regime 
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In the impulse approximation regime the nuclear response to a weakly interacting probe can be written in terms 
of the measured nucleon structure fuctions and the target spectral function, yielding the energy and momentum 
distribution of the constituent nucleons. We discuss a calculation of charged current neutrino-oxygen interactions 
in the quasielastic channel, carried out within nuclear many body theory. The proposed approach, extensively and 
successfully employed in the analysys of electron-nucleus scattering data, allows for a parameter free prediction 
of the neutrino-nucleus cross section, whose quantitative understanding will be critical to the analysis of the next 
genaration of high precision neutrino oscillation experiments. 



1. Introduction 

In experimental searches of oscillations, neutri- 
nos are detected trough their interactions with 
nucleons bound in nuclei (e.g. oxygen). There- 
fore, the knowledge of the neutrino-nucleus cross 
section at a fully quantitative level will be critical 
to reduce the systematic uncertainty of the next 
generation of high precision measurements. 

Accurate calculations of the weak nuclear re- 
sponse at moderate momentum transfer can be 
carried out within nuclear many body theory 
(NMBT) using nonrelativistic wave functions to 
describe the target initial and final states and 
expanding the current operator in powers of 
(|q|/m), q and m being the momentum transfer 
and the nucleon mass, respectively. Applications 
of this approach to the analysis of electron scat- 
tering data has been very successful (for a review 
see, e.g., ref.pP). 

At the higher values of |q| corresponding to 
energies of the beam particles in the few GeV re- 
gion, the description of the final states in terms of 
nonrelativistic nucleons is no longer possible, and 
some simplifying assumptions have to be made in 
order to take into account both the presence of 
relativistic particles and the occurrence of inelas- 
tic processes. 



An approximation scheme widely employed to 
describe the region of high momentum transfer 
is based on the idea that, as the beam particle 
probes a region of extension ~ l/|q| of the target, 
at large |q| the scattering process involves only 
one nucleon, the remainig A-l particles acting as 
spectators. 

In the simplest implementation of this scheme, 
usually referred to as impulse approximation 
(IA), the hadronic state produced at the weak 
interaction vertex is assumed to be totally decou- 
pled from the recoiling system, so that the de- 
scription of its motion reduces to a purely kine- 
matical problem. On the other hand, the target 
intial state, as well as the final state of the specta- 
tor particles, can be safely treated within NMBT 
using nonrelativistic dynamics. 

In Section 2, after giving the expression of the 
neutrino-nucleus cross section in the IA regime, 
we briefly outline a theoretical approcah, based 
on NMBT, that allows for both a realistic de- 
scription of the initial state and the inclusion of 
corrections arising from final state interactions 
(FSI). The results of calculation for the test case 
of quasielastic charged current interactions with 
oxygen are presented in Section 3, while Section 
4 is devoted to a summary and the conclusions. 
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2. The neutrino-nucleus cross section in 
the impulse approximation regime 

The Born approximation cross section of the 
weak charged current process 



vi + A 



X 



(1) 



where A and X denote the target nucleus and the 
undetected hadronic final state, respectively, can 
be written in the form 
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In the above equation, G is the Fermi coupling 
constant, Eg is the energy of the final state lep- 
ton and k and k' are the neutrino and charged 
lepton momentum, respectively. The tensor L M „ 
is fully specified by the kinematical variables of 
the leptons, while W^ v is defined in terms of the 
target initial and final states and current accord- 
ing to 

= J2( A \ J ^\ X )( X \ JV \ A ) 

xS^(P A + q-P x ) , (3) 



X 



Pa, Px and q being the target initial and final 
four momenta and the four momentum transfer, 
respectively. 

The IA scheme outlined in Section 1 is based 
on the assumptions that: i) the target nucleus is 
seen by the neutrino as a collection of individual 
nuclcons and ii) the hadron produced at the weak 
interaction vertex is unaffected by either Pauli 
blocking or final state interactions (FSI) with the 
spectator nucleons. 

According to this picture, the state \X) appear- 
ing in Eq.Q can be written in a simple factor- 
ized form and the tensor W^ v of Eq.© reduces 
to (see, e.g. ref.[3]) 



d 4 pP(p) uf(p t q) , 



(4) 



where P(p) (p = (po,p)) is the nuclear spectral 
function [2], yielding the probability of removing 
a nucleon of momentum p from the target leav- 
ing the residual system with energy E = m — po. 
Eq.Q shows that the IA formalism allows one 



to describe scattering off a bound nucleon carry- 
ing four momentum p in terms of the tensor w^ v , 
describing the weak charged current interactions 
of a nucleon in free space. Binding effects are 
taken care of by replacing the bound nucleon four 
momentum p = (po, p) with p = (-\/p 2 + m2 > P) 
and the four momentum tranfer q = (v, q) with 
q=p-p + q. [51315] . 

The IA cross section can be obtained from 
Eqs.J2J and 10} using the available parametriza- 
tions of the structure functions (both clastic and 
inelastic) entering the definition of w^ v and the 
nuclear spectral functions resulting from calcula- 
tions carried out within nuclear many body the- 
ory (NMBT) 0. 

NMBT also provides a consistent framework to 
include corrections to the IA in a systematic fash- 
ion. The effects of FSI, that have long been rec- 
ognized to be large in both inclusive and exclusive 
electron-nucleus scattering processes (a series of 
review papers on electron-nucleus scattering can 
be found in ref.pQ), can be taken into account 
within the approach of ref. [3| , based on the as- 
sumptions that i) the hadron produced at the 
weak interaction vertex moves along a straight 
trajectory with constant velocity (eikonal approx- 
imation) and ii) the spectator nucleons are seen 
as a collection of fixed scattering centers (frozen 
spectators approximation). The resulting cross 
section can be written in terms of the IA results 
according to 



.(5) 
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d^h = J dE ' e h{Ei ~ K) (^S^ 

The above equation shows that the the occurrence 
of FSI leads to a redistribution of the strength dic- 
tated by the folding function f q (E), whose shape 
is strongly affected by nucleon-nucleon correla- 
tions 0. In absence of FSI f q (E) -> 5(E) and 
the IA cross section is recovered. 

3. Results 

We have computed the cross section of the 
quasielastic process 

v e + w O-^e- +p + X , (6) 

where X denotes the recoiling spectator system. 
The spectral function apperaing in Eq.(0J has 
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been obtained within the Local Density Approx- 
imation (LDA) [S], while for the Dirac, Pauli 
and axial form factors entering the definition of 
the nucleon tensor w^ v we have used the simple 
dipole parametrization. 

The folding function has been evaluated using a 
parametrization of the measured nucleon-nucleon 
scattering amplitude [i] and Monte Carlo configu- 
rations sampled from the probability distribution 
associated with the realistic oxygen wave function 
ofref. H]. 
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Figure 1. Differential cross section da /d£l e dE e for 
neutrino energy E = 1 GeV and electron scattering 
angle 8 e = 30° . The IA results are represented by the 
dashed line, while the solid line correspond to the full 
calculation, including the effects of FSI. The dotted 
line shows the prediction of the FG model with Fermi 
momentum k F = 221 MeV and average separation 
energy e = 35 MeV. 



Fig. ^shows the cross section corresponding to 
neutrino energy E = 1 GeV, and electron scat- 
tering angle 9 e = 30°, plotted as a function of 
the energy transfer E — E e . Comparison between 
the solid and dashed lines shows that the inclu- 
sion of FSI results in a sizable redistribution of 
the IA strength, leading to a quenching of the 
quasielastic peak and to the enhancement of the 
tails. For reference, we also show the cross sec- 
tion predicted by the Fermi gas (FG) model with 
Fermi momentum kp — 221 MeV and average 
separation energy e = 35 MeV. It appears that 



nuclear dynamics, neglected in the oversimplified 
picture in terms of noninteracting nucleons, plays 
a relevant role. 

It has to be pointed out that the approach of 
ref . [3] , while including dynamical correlations in 
the final state, does not take into account statis- 
tical correlations, leading to Pauli blocking of the 
phase-space available to the outgoing nucleon. 

A rather crude prescription to include the ef- 
fect of Pauli blocking amounts to modifying the 
spectral function according to 

P(p)^P(p)6(\p + q\-k F ) (7) 

where k F is the average nuclear Fermi momen- 
tum, defined as 

k F = J d 3 r PA {v)k F (r), (8) 

with fc F (r) = (3vr 2 / 9A(r)/2) 1 / 3 , p A (r) being the 
measured nuclear density distribution. Eq.lpjjl 
yields k F = 209 MeV. Note that, unlike the spec- 
tral function, the quantity defined in Eq.Q docs 
not describe intrinsic properties of the target only, 
as it depends explicitely upon the momentum 
transfer. 
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Figure 2. Differential cross section da/dQ 2 for neu- 
trino energy E = 1 GeV. The dot-dash line shows 
the IA results, while the solid and dashed lines have 
been obtained using the modified spectral function 
ef Eq.lQ, with and without inclusion of FSI, respec- 
tively. 
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The effect of Pauli blocking is hardly visible in 
the differential cross sections shown in Fig. as 
the kinematical setup corresponds to Q 2 ~ 2.4 
GeV 2 at the quasielastic peak. On the other 
hand, it becomes very large at lower Q 2 . 

Fig. [3 shows the calculated differential cross 
section da/dQ 2 for neutrino energy E = 1 GeV. 
The dot-dash and dashed lines correspond to the 
IA results with and without inclusion of Pauli 
blocking, respectively. It clearly appears that the 
effect of Fermi statistic in suppressing scattering 
at Q 2 < .2 GeV 2 is very large and must be taken 
into account. The results of the full calculation, 
in which FSI are also included, are displayed as a 
full line. The results of Fig. [21 suggest that Pauli 
blocking and FSI may explain the deficit of the 
measured cross section at low Q 2 with respect to 
the predictions of Monte Carlo simulations [H]. 

The energy spectrum of the produced electrons 
obtained using the modified spectral function of 
Eq.0 is displayed in Fig. The solid and 
dashed lines corresponds to the results of calcu- 
lations carried out with and without FSI. As in 
Fig. ^ it appears that inclusion of FSI results 
in a redistribution of the strength, leading to a 
quenching in the region of the peak and to the 
appearance of a tail at large energy, extending 
all the way to the boundary of the kinematically 
allowed region. 

4. Conclusions 

NMBT provides a parameter free and computa- 
tionally viable framework to study the nuclear re- 
sponse to weakly interacting probes, whose quan- 
titative understanding will be critical to the anal- 
ysis of the next genaration of high precision neu- 
trino oscillation experiments. 

In the IA scheme, which is likely to be ap- 
plicable to describe scattering of few GeV neu- 
trinos, NMBT allows for a consistent treatment 
of dynamical nucleon-nucleon correlations, that 
largely affect both the nucleon spectral function 
and FSI. Electron scattering experiments have 
provided overwhelming evidence that the inde- 
pendent particle model fails to explain the data 
and inclusion of correlation effects, leading to a 
sizable redistribution of the strength, is indeed 
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Figure 3. Energy spectrum of the produced elec- 
trons in charged current interactions at neutrino en- 
ergy E — 1 GeV. The solid and dashed lines have 
been obtained including the effect of Pauli blocking 
with and without FSI, respectively. 



required. 

As a final remark, it has to be emphasized 
that the extension of the approach described in 
this paper to include exclusive as well as inelas- 
tic channels, along the lines of the work done for 
electron-nucleus scattering, is straightforward. 
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